Glass breaks like metals, but at the nanometer scale 
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We report in situ Atomic Force Microscopy experiments which reveal the presence of nanoscale 
damage cavities ahead of a stress-corrosion crack tip in glass. Their presence might explain the 
departure from linear elasticity observed in the vicinity of a crack tip in glass. Such a ductile fracture 
mechanism, widely observed in the case of metallic materials at the micrometer scale, might be also 
at the origin of the striking similarity of the morphologies of fracture surfaces of glass and metallic 
alloys at different length scales. 

PACS numbers: 62.20.Mk, 81.40.Np, 87.64.Dz 



Glasses are the most common example of "brittle" ma- 
terials which break abruptly, without first deforming in 
an irreversible way as metals do. In metallic alloys, cracks 
usually progress through the coalescence of damage cav- 
ities, which nucleate within microstrutural defects (sec- 
ond phase precipitates, grain boundaries...) or at the in- 
terface between the matrix and the heterogeneities. This 
"ductile" fracture mode,widely observed for a large va- 
riety of metallic alloys p|, leads to very rough fracture 
surfaces which have been extensively studied over the 
last eighteen years 0, 0|- F° r s l° w crack propagation, 
the fracture surfaces of glass may appear very flat (in 
the so-called "mirror" zone) if examined with an optical 
microscope. However, when analysed at the nanometer 
scale with an Atomic Force Microscope (AFM), they re- 
veal a roughness which is strikingly similar to the one 
exhibited by metallic fracture surfaces 0, 0|- The only 
difference actually resides in the length scales involved, 
which are several orders of magnitude smaller in the case 
of glass. Could it be because, despite the conventional 
belief, the fracture mechanisms of these two categories of 
materials are similar, although taking place at different 
length scales? We report here the first experimental ev- 
idence of a ductile fracture mode in a vitreous material 
at a temperature much lower than the glass transition 
temperature Tg. Such scenario was predicted by Molec- 
ular Dynamics simulations 0, 0|, but it had never been 
observed experimentally up to now. The experiments re- 
ported here clearly show that slow fracture in glass pro- 
gresses through the nucleation, growth and coalescence 
of damage cavities at the nanometer scale. These cavi- 
ties are shown to be correlated to the non-linear elastic 
zone observed in the vicinity of the crack tip Q. Possible 
origins of the cavities' nucleation are conjectured and the 
consequences of such ductile fracture mode in glass are 
discussed. 

Experimental setup. The experimental setup is illus- 
trated in Fig. ^ All the experiments are performed at a 
constant temperature of 22.0 ±0.5 °C and in a leak-proof 
chamber under an atmosphere composed of pure nitro- 
gen and water vapor at a relative humidity of 42 ± 1% 



after preliminary out-gassing. Fracture is performed on 
DCDC Q (Double Cleavage Drilled Compression) paral- 
lelepipedic (4 x 4 x 40 mm 3 ) samples of aluminosilicate 
glass. A thermal treatment (660°C) was performed be- 
fore each fracture experiment in order to remove residual 
stresses Q. The 4 x 40 mm 2 surfaces are optically pol- 
ished (RMS roughness is 0.25 nm for a 10 x 10 /im 2 scan 
size). In the centre of two parallel 4 x 40 mm 2 surfaces 
and perpendicularly to them a cylindrical hole (radius 
a = 0.5 mm) is drilled. Its axis defines the z-direction. 
The x-axis (and y-axis) are parallel to the 40 mm (and 
4 mm) long side of the 4 x 40 mm 2 surface. A compressive 
load is applied perpendicularly to the 4x4 mm 2 surfaces. 
The external stress a is gradually increased by the slow 
constant displacement (0.02 mm/min) of the jaws of the 
compression machine. Once the two cracks (symmetric 
to the hole axis) are initiated, the jaws 'displacement is 
stopped. The crack front then propagates along the x- 
axis in the symmetry plane of the sample parallel to the 
(x,z) plane. In this geometry, the stress intensity factor 
Ki is given ,S| by: Kj = a^/a/ (0.375c/ a + 2), where c is 
the length of the crack (Fig.QJi). 

At the very first moments, the crack propagates very 
quickly. In this regime, the crack velocity v is independent 
of the chemical composition of the surrounding environ- 
ment (l0( |. As the crack length c increases, Ki decreases, 
and v decreases quickly. Under vacuum, the crack stops 
for Ki smaller than a critical value Kj c referred to as the 
toughness of the material. But in a humid atmosphere, 
the corrosive action of water on glass allows for slow 
crack propagation at much lower values of the stress [Icj . 
The crack motion within the external (x,y) sample sur- 
face is then slow enough to be monitored by our exper- 
imental system combining optical microscopy and AFM 
(Fig. Qb). Optical image processing gives the position of 
the crack tip and consequently the "instantaneous" ve- 
locity for v ranging from 10 -6 to 10~ 9 m.s -1 . By AFM 
measurements - performed in a high amplitude resonant 
mode ("tapping" mode) -, one probes the crack tip neigh- 
bourhood at magn ifications ranging from 75 x 75 nm 2 
to 5 x 5 /im 2 1 1 1 and the crack tip motion at veloci- 
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ties ranging from 10 to 10 12 m.s . The data pre- 
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sented below are obtained for Ki = 0.43MPa.m 1/z and 
v = 3.10- 11 m/s. 

Evidence of nano-scale damage cavities. Typical topo- 
graphical frames in the neighbourhood of the crack tip 
are presented in Fig. [21 They clearly reveal cavities of 
typically 20 nm in length and 5 nm in width ahead of 
the crack tip (Fig. EK)- These cavities grow with time 
(Fig. EJd) until they coalesce (Fig. Eh). 

To ensure that the spots observed ahead of the crack 
tip are actually damage cavities which grow further and 
coalesce with the main crack leading to failure, we use the 
Fracture Surface Topography Analysis (FRASTA) tech- 
nique [12] first introduced to study damage in metallic 
alloys. It consists in analyzing the mismatch between the 
two fracture surfaces. To understand how the FRASTA 
technique can provide physical details on the fracture 
mechanism, let us now consider how a ductile material 
breaks: In such a medium, the load application first in- 
duces a local plastic flow before generating any local fail- 
ure at the level of the stress concentrators. Then, when 
local failure actually occurs, the stress applied on the 
newly formed void surfaces vanish, the applied load is 
redistributed to nearby unbroken material, and plastic 
deformation is not undergone anymore at the level of 
these void surfaces. Consequently, in such a ductile sce- 
nario, each cavity initiation is accompanied by local irre- 
versible plastic deformations printed in relief on the devel- 
oping fracture surfaces (the crack lines when the method 
is applied in two dimensions as in the present case) that 
should remain visible after the cavities have coalesced 
and the crack has crossed. 

The method consists in placing the upper fracture sur- 
face (crack line in 2D) under the lower one until no void 
is left. Then, the two surfaces (or lines) are pulled away 
from each other along the direction perpendicular to the 
fracture plane (or direction). This is what happens dur- 
ing the fracture process when the external strength is 
applied at a constant displacement rate. At the small 
length scales considered here, this assumption can be 
made. Cavities therefore appear in the chronological or- 
der. 

The crack lines are first determined by binarising the 
image of the sample after fracture (Pig.Et)j an d the un- 
broken material is reconstituted virtually by placing nu- 
merically the lower line over the upper one (Fig.^b). The 
lower crack line is gradually translated in the direction of 
decreasing y (Fig. |2h>) , and the cavities appear and grow 
in the chronological order. The structure obtained for a 
given displacement, i.e. at a given time, is superimposed 
on images recorded prior to failure and shown to cor- 
respond actually to cavities observed at this given time 
(Fig.Et). This indicates that the spots observed prior to 
failure are indeed depressions which are marked in relief 
on the final crack, and hence, are actual damage cavities. 

Displacement field. A consequence of this "nano-scale 
ductility" can be seen in the displacement field around 
the crack tip (see also references 0] for related discus- 



sion): For a slit-like plane crack in an ideal Hookean 
continuum solid, the stress components cr^ at a given 
point M whose cylindrical polar coordinates are (r, 9, z) 
in the vicinity of the crack tip take the form Oij (r, 9, z) — 
fij(9)Ki^/27rr where the functions are completly de- 
termined from linear elastic theory 13] . Moreover, the 
z-displacement u z of the specimen surface (flat when no 
load is applied) is given by u z (r, 6) = f. e zz (r, 6, z)dz 
where h = 4 mm is the thickness of the specimen and 
e zz is the direct strain in the z-direction. If the linear 
elastic stress-strain relation held on the nanometer scale, 
one gets e zz = (a zz - v(a xx + <J yy ))/E (where E and 
v are respectively the Young Modulus and the Poisson 
coefficient of the material) and consequently u z (r, 8) = 
G(0)hvKi I (E\J 2nr) , where G{9) is a function equal to 
G{9) = 2cos(6»/2) for mode I fractures Q. 

Measurements of u z profiles have been performed on 
lxl fim 2 AFM topographical frames (Fig. 0k) along the 
direction of crack propagation (Fig. QJd) and perpendic- 
ularly to it (Fig. For both profiles, u z departs from 
the linear elastic r -1 / 2 scaling for r smaller than a thresh- 
old r c highly dependent on 9: for 9 = 0°, r c — 100 nm 
while for 9 = 90°, r c — 20 nm. These short-range depar- 
tures from the linear elastic behaviour may be related 
to the presence of cavities although many other phe- 
nomena could be responsible for this discrepancy. How- 
ever, the fact that the order of magnitude of the ratio 
r c (9 = 0°)/r c (9 = 90°) - much higher than the ratio of 
the cosine terms in the linear elastic expression of u z - is 
close to the aspect ratio of the observed damage cavities, 
strongly suggests a correlation between damage and non 
linear elasticity. 

Discussion. Similar fracture experiments performed in 
amorphous Silica specimens reveal similar damage cavi- 
ties. This suggests that their existence does not depend 
on the precise chemical composition of the studied glass. 
The cavities' nucleation should be found more likely in 
the amorphous structure, which contains inherent atomic 
density fluctuations at the nanometer scale. Such atomic 
density fluctuations have been evidenced by van Brutzel, 
Rountree, et al. |j in the structure of simulated amor- 
phous silica by Molecular Dynamics: The Si and O atoms 
are shown to form Silica tetrahedra connected together to 
build rings of different sizes ranging from 3 to 9 tetrahe- 
dra. At larger length scales, ranging from 1.5 nm to 6 nm, 
the density of these rings is found to fluctuate with high 
density areas surrounded by low density areas. Moreover, 
the Molecular Dynamics of van Brutzel Q show that, at 
this length scale, crack propagates by growth and coa- 
lescence of small cavities which appear in areas with low 
density of rings, ahead of the crack tip. They behave as 
stress concentrators and grow under the stress imposed 
by the presence of the main crack to give birth to the 
cavities actually observed in the AFM frames. 

Another scenario was also proposed by Falk [TH Hflj: 
A small change in the interactomic potential was shown 
to be sufficient to generate a ductile to brittle transi- 
tion in a 2D polydisperse packing of beads interacting 
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through Lennard- Jones potential ^^|. This sensitivity 
was related to the existence of non-affme rearrangements 
zones (groups of atoms whose motion cannot be approx- 
imated by a linear local strain tensor) whose density de- 
pends strongly on the interatomic potential range. Un- 
fortunalty, the AFM resolution is not presently sufficient 
to directly evidence these non-affme molecular rearrange- 
ments. 

Moreover, two roughness regimes have been actually 
observed on post-mortem glass fracture surfaces 0, E|- 
The low length scales regime (from 1 nm to a crossover 
length £ c ranging from 10 to 30 nm depending on the av- 
erage crack velocity) was interpreted as the consequence 
of the extension of an isolated damage cavity; the regime 
at larger length scales (from £ c to 100 nm) is thought to 
be linked to the structure composed of the various cor- 
related damage cavities fljj • The order of magnitude of 
£ c is in good agreement with our present observations, 
where cavities at coalescence are a few tens of nanome- 
ters wide. 

Here, let us note that AFM observations are performed 
on the sample surface, where the mechanical state is dif- 
ferent from that of the bulk. It must be remembered that 
the stress components at the surface are contained within 
the free surface plane, while the deformations are fully 
three-dimensional. On the contrary, the bulk is in a con- 
dition of plane strain and fully three-dimensional stresses. 
Hence, the observed sizes and growth rates of cavities at 
the surface may well differ from those in the bulk. New 
experiments using the FRASTA method in three dimen- 
sions applied to the post-mortem study of the fracture 
surfaces are currently being performed, in order to have 
access to the three dimensional structure of bulk damage 



and its evolution. Through this new set of experiments, 
one should be able to correlate also the damage structure 
to the fracture surface morphology. 

Moreover, the structure of damage, which influences 
macroscopic mechanical properties such as fracture 
toughness and lifetime, should then be linked to the glass 
composition and nanostructure [fJ,ElE|]- Complemen- 
tary analyses addressing the question of the chemical 
bonds on the fracture surface will also be performed. 

Finally, the fact that glass at temperatures far below 
the glass transition temperature Tg, considered up to now 
as the archetype of pure brittle material, joins the class 
of damageable materials should have important conse- 
quences for its mechanical properties. In applications, 
the design of structures using glass might be modified 
to take this behaviour into account, especially for slow 
crack propagation processes. 

The similarity between the damage modes of materials 
as different as glass and metallic alloys is an important 
clue to understand the origin of puzzling universal be- 
haviors, hence sheding new light into the basic physical 
mechanisms of fracture. 
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FIG. 1: experimental setup, (a): Sketch of the DCDC geometry (b): Picture of the experiment. 
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PIG. 2: Sequence of successive topographic AFM frames showing the crack propagation at the surface of the specimen. The 
scan-size is 75 x 75 nm 2 and the heights range over 2 nm. The recording time for one frame is around 3 mn and two successive 
frames are separated by 30 mn. The crack front propagates from the left to the right (^-direction) with a mean velocity v of 
around 10 -11 m/s. (a): evidence of nanometric damage cavities before the fracture advance, (b): growth of the cavities, (c): the 
crack is advancing via the coalescence of all the cavities 



FIG. 3: Fracture Surface Topographical Analysis (FRASTA). (a) Frame 0J; (broken sample) is binarised and the contours of 
the crack are determined, (b) The lower line is first numerically raised over the upper one and then gradually displaced in the 
direction of decreasing y, as schematized by the arrow. Cavities are coloured in red. (c): Result of the method: superimposition 
of the obtained cavities on the image recorded prior to complete failure. 



FIG. 4: Measurements of the surface deformations and comparison with the predictions for an ideal Hookean material. The crack 
propagates from left to right (x positive), (a): typical AFM topographical frame of the vicinity of the crack tip. The scan-size 
is 1 x 1 fim 2 and the heights range over 3 nm. The white vertical (respectively horizontal) dotted line sets the ^-coordinate 
Xf (respectively the y-coordinate yf) of the crack tip. (b) (respectively (c)) top: Plot of the z profile along (respectively 
perpendicularly to) the direction of crack propagation. The open circles correspond to experimental data while the full line 
corresponds to the prediction z = Zo — Ar~ 1 ^ 2 (where zo and A are fit parameters) given for an ideal Hookean solid. Bottom: 
Log-log plot of the displacement u z = zo — z versus the distance r = x — Xf (respectively r = y — yf) from the crack tip. For 
r < r c , the 5-profile departs from the predictions given by the linear elastic theory. 
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